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We study the local density of states (LDOS) of two-dimensional electrons in the presence of 
spin-orbit (SO) coupling. Although SO coupling has no effect on the average density of states, 
it manifests itself in the correlations of the LDOS. Namely, the correlation function acquires two 
satellites centered at energy difference equal to the SO splitting, 2lo S o, of the electron Fermi sur- 
face. For a smooth disorder the satellites are well separated from the main peak. Weak Zeeman 
splitting luz "C wso in a parallel magnetic field causes an anomaly in the shape of the satellites. We 
consider the effect of SO-induced satellites in the LDOS correlations on the shape of the correlation 
function of resonant-tunneling conductances at different source-drain biases, which can be measured 
experimentally. This shape is strongly sensitive to the relation between u>so and u> z . 

PACS numbers: 73.40.Gk, 73.20.Fz, 71.70.Ej 



The density of states, is , in 2D electron gas is energy-independent. It remains energy-independent 
in a parallel magnetic field, which causes the spin-splitting, 2uj z , of the electron spectrum. The 
density of states, corresponding to each spin branch, is equal to v^jl. Spin-orbit (SO) coupling also 
results in the splitting of the electron spectrum into two branches [1]. In this case the density of 
states in each branch, fj, — ±1, depends on energy 



where ±co> so is the splitting of a state which had the energy e ^> u so in the absence of the SO- 
coupling. As can be seen from Eq. (1), the net density of states is still identically equal to u . Since 
the magnitude of the splitting, 2u so , at the Fermi level e = E F is inversely proportional to the length 
of the spin rotation, this magnitude represents an important characteristics of the 2D structure, 
containing the electron gas. The importance of the spin- rotation length was appreciated since 1990, 
when the proposal for device application of the spin-polarized currents was put forward [2]. This 
proposal has lately attracted a lot of interest. 

If the magnitude of the SO-splitting is large enough, it can be inferred from the beating pattern 
of the Shubnikov-de Haas oscillations, as was first demonstrated by Dorozhkin and Olshanetskii [3] 
for 5*i-based structures, and subsequently [4] by Luo et al. for narrow-gap heterostructures. While 
for large u; so (several meV), as in [4], the beats of the Shubnikov-de Haas oscillations yield a rather 
accurate value of SO-splitting, extracting small splittings from the beating pattern is complicated in 
two regards: (i) very low magnetic fields (with cyclotron energy smaller than u so ) and, correspond- 
ingly, very low temperatures are required to observe the beatings; (ii) presence of even moderate 
disorder suppresses the Shubnikov-de Haas oscillations by the Dingle factor, which gets large at low 
fields. 

In the present paper we demonstrate that the disorder can, actually, reveal the SO-splitting even if 
u; so is smaller than the single particle scattering rate, r~ x . In a disordered sample the local density 
of states (LDOS) fluctuates randomly in space and these fluctuations contain information about SO 
coupling. Our main point is that this information, which is lost in the average LDOS, is preserved in 
the correlation function of LDOS at two different energies, P(ei, e^)- This function, whose magnitude 
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is inversely proportional to the disorder, exhibits peaks at (ei — 62) = ±2u so . Even for uj so t <C 1 the 
peaks are well pronounced when the disorder is smooth, so that the transport relaxation time, r tr , 
is much longer than r. This is illustrated in Fig. 1. Smooth disorder implies that the momentum 
transfer in a single scattering act is small compared to the Fermi momentum, k F . In this way, the 
states with energies e + u so and e — u so remain correlated after many, ~ T tr /r, scattering acts. 
Therefore, the SO-induced satellites in P(e\, £2) are sharp even in the presence of a strong disorder, 
t -1 > u so , if the condition u so r tr ^> 1 is met. In terms of underlying physics, the observation that a 
small SO-splitting can be resolved in a smooth potential is quite analogous to the observation made 
in Ref. [5] concerning the manifestation of the Landau levels in P(e\, €2) in the presence of a strong 
disorder, when the Landau level structure in the average density of states is completely smeared out 
[6]. Let us emphasize, though, that in our case (unlike Ref. [5]) the density of states in the absence 
of disorder is constant, which does not contain u so . Thus, it is only in the presence of disorder that 
LDOS becomes sensitive to the SO splitting. 

Concerning the experimental consequences of SO-induced satellites in P(ei,€2), we note that the 
local density of states governs the probability of electron tunneling from the impurity in the barrier 
into the 2D gas [7-9]. This fact was recently utilized in the resonant-tunneling spectroscopy [10-14] 
of the LDOS. More specifically, the correlator of the LDOS determines the behavior of correlation 
function, C(5V), of the fluctuations of the tunneling conductance with the change, 8V, of the 
source-drain bias. From the analysis of C(SV), measured experimentally, the authors [13] were able 
to extract the quantitative information about the structure of states and their inelastic lifetime in 
the disordered emitter. In the present paper we demonstrate that the SO-induced correlations in 
the LDOS give rise to the new characteristic features in C(5V). We also demonstrate that these 
new features are extremely sensitive to a weak parallel magnetic field. 

The paper is organized as follows. In Sect. II we derive the analytic expression for P{e\ — 62) in the 
case of a smooth disorder. In Sect. Ill we give the qulitative explanation for the ratio between the 
amplitudes of the main peak in P(t\ — e 2 ) and of the SO-induced satellites. In Sect. IV we use the 
shape of P{e\ — £2) to calculate the SO-induced features in the correlator, C(5V), of the tunneling 
conductances. Concluding remarks are presented in Sect. V. 



II. CALCULATION OF P(ei - e 2 ) 
A. Hamiltonian and eigenfunctions 

We choose the conventional form, a(cr x k)n , of the SO-term originating from the confinement 
potential asymmetry [1]. Here a is the coupling constant, cr is the spin operator, and n is the unit 
vector normal to the two-dimensional plane. In a parallel magnetic field, that induces the Zeeman 
splitting 2u z , the Hamiltonian of a free electron is 

h 2 k 2 

H = h a(& x k)n + uo z a X) (2) 

2m 

where m is the electron mass. The eigenfunctions of the Hamiltonian Eq. (2) are classified according 
to the chirality, \x = ±1, and form two branches of the spectrum, as illustrated in Fig. 1. In the 
vicinity of the Fermi surface the spectrum, E^(k), can be simplified 

E li (k) = E F + e li {k), (3) 

where e /i (fe) is defined as 
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e^k) = hv F (k- k F ) + fitt(k), (4) 

and 



Q(k) = y/u% + Jl + 2u so uj z sin0 fe . (5) 

Here v F = hk F /m is the Fermi velocity, u so = ak F , and <pu is the azimuthal angle of k. It is convenient 
to introduce the projection operators A M (fe), which are defined through the eigenfunctions 



X ^ 2 1 / 2 ( -iiiexp(i(p h ) ) ' 



(6) 



of the Hamiltonian Eq. (2) as A /1 (fe) = so that [15] 



A„(fe) = £( • \. ^«P(-W|, (7) 

where the angle fk is related to the azimuthal angle 4>k as 

tan tp k = tan fc H — — — . (8) 

U so COS fc 

In terms of projection operators the Hamiltonian can be presented in a simple form H = 
^E,(k)A„(k). 

B. Correlator of the LDOS in the chirality representation 

In the presence of disorder, the LDOS at energy e at point r is defined as 

u(r, e) = ^-Tr [d A (r, r, e) - G R (r, r, e)l , (9) 

where the advanced and retarded Green's functions are 2x2 matrices in the chirality space, and the 
trace is taken over different chiralities. The main contribution to the LDOS correlator comes from 
the cross terms G^G^ 



P(e-e') = ^(5v(r,e)5v(r,e>)) = -^R e (TrG>, r, e)TrG A (r, r, e')>, (10) 



where 5u(r, e) is the deviation of the LDOS from the disorder-free value, v Q . It is easy to see from 
Eq. (4) that vq = m/irh 2 , and does not depend on energy. To evaluate the correlator Eq. (10) it is 
convenient to use the chirality representation of the Green functions 

where r is the scattering time. As in the absence of the SO-coupling, the LDOS correlator rep- 
resents the sum of the diffuson and cooperon contributions [16]. Generalized expressions for the 
corresponding contributions in the chirality representation read 

x G£ (e, p + q)Gl (e, p> + q)G A 2 (e\ p)G A (e', p') , (12) 
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xG«(e,q-p)G*(e,q-p')G^( e ',p)G*(e',p')= P D (e' -e), (13) 

where 

^^,(p,p') = xt(p)x v (p') = ^{l + /^exp(i[v? p - w])}. (14) 

The two-particle vertex functions r^^(p,p', q, e) in Eqs. (12), (13) satisfy the matrix Dyson-type 
equation 



I™(p,p',<Z,e) = S(p-p')<S>^ Vl (p,p')% 2 „ 2 (p,p') £ K^ll{p^q^)Y^ 2 {p u p', q ,e) 

(15) 

with a kernel 

*£H(P,Pi, 9, 6) = $„ 1>€l (p,p')^ >6 (p,l/)G ! g(a; + c,pi + <z)G f A >, Pl )S(p - pf), (16) 

where the function S(p — p') is the Fourier transform of the correlator of the random potential. 
Upon integration over \pi\ Eq. (15) takes the form 

r»P,p',<?,e) = s(p- j/)<iw(p>jO<«W(p>P') ^ ^S§(P.Pi.9.e)rgS(Pi,P',9,c) 

(17) 

where we have introduced the modified kernel , defined as 

K^(v v a e) = (—) ®^,UPiP')®^2(PiP')S{p-p') m 

w& V 2tt 7 1 - i( w - eifi(pi) + 6«(Pi + + ^9«pT cos(0 p - q ) ' 1 ^ 

When the random potential is smooth, the function S restricts the difference p—p' within a narrow 
domain \p — p'\ <C k F . As we will see below this leads to a drastic simplification of the system Eqs. 
(17), (18). 

C. Smooth potential; |ei — £2! ~ ^so 

In the case of a smooth potential the factors $ M) ^(p, p') for coinsiding chiralities, /x = 77, differ 
strongly from those with opposite chiralities, /x = —77. Indeed, when |p — p'| is small compared to 
k F , we have P ~ P ' and thus y?p ~ yy, so that <C 1, whereas is close to unity and can 

be presented as 

^Jp,p') = l-^^. (19) 

Using Eq. (8), the factor $ MiAt can be expressed through the angle, 4> p — p /, between the vectors p 
and p' as follows 

%AP,P') = 1 " " (20) 
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The latter simplification allows to set &n,-n = and = 1 everywhere except for the kernel of 
the Dyson equation. As a result, the system Eqs. (17) gets decoupled into closed equations for the 
elements, with a ^ rt, of the matrix r^!£. The underlying reason for this decoupling is that 
these elements describe the two-particle motion, in course of which one particle moves within the 
branch \x whereas another particle moves within the branch 77. For a smooth potential coupling of 
these elements to the other branch-nonconserving elements is small in parameter (1 — <C 1. 

With the above simplification Eq. (17) takes the form [17] 



r%M,q,u) = S{K\4> P -<t>p\) +/^^(p,Pi,9^)r^(pi,p',q,a;), (21) 

while the kernel Eq. (18) simplifies to 

K™(v p, Q uj)=(—\ l^,>,Pi)| 2 £(M0 P -0 P J) m) 

™ lP ' PU 9 ' ' V 2vr A-ir - v) «(p) - cos(0 p - g )] ' 1 J 

Substituting Eq. (22) into Eq. (21) yields an integral equation for the vertex T^. Note, that this 
equation differs from the standard equation for the vertex in the absence of the SO-coupling only by 
the factor |<E» Mj/1 (p, pi)\ 2 in the integrand, where $ MiM (p, p\) is defined by Eq. (19). Thus, if the small 
difference of this factor from unity is neglected, the solution for would contain a conventional 
diffusive pole at u — (/j, — rf)Q(p) — iDq 2 , where D = v 2 T tr /2 is the diffusion coefficient, and 

V 1 = \j-^U>>- S (\ k ~ Pi) W*0 - E,{p)) . (23) 

is the transport relaxation time. Taking into account the small correction originating from the 
difference 1 — Pi)| 2 amounts to the imaginary shift of the pole position by ir~J, where r int 

is defined as 

r-l(k) = ll^y 2 (^- ^? S(\k - p\)S(E,(k) - E,(p)) . (24) 

The meaning of Ti nt is the interbranch scattering time. Indeed, a general expression for the scattering 
time between the branches with different chiralities can be written as 



^U k ) = / (Mfc)A-») S(\k-p\)5(E,(k) - E_,(p)) 



I 



dp 

2^ 
dp 

(2^F 



(l - |$,, M (fc,p)| 2 ) S{\k - p\)5(E,(k) - E_,(p)) ■ (25) 



Comparing Eqs. (24) and (25), we see that they differ only by the arguments of the (^-functions. 
This difference is negligible when u so <C E F , and thus we have T^-^k) m Ti nt (k). We also see that 
the integral Eq. (24) in the expression for differs from the integral in Eq. (23) for r^ 1 only by 
the replacement — > Using the relation Eq. (20) between the two angles, we can express r int 
through the transport relaxation time 

r,M) = 2Ttr M±^±^i^l!. (26) 

^so (^so +w z sm0) 

With shifted diffusion pole, the final expression for the vertex takes the form 

T%M,q,»)= , , ^yrfe 1 ?. . - .r (27) 



u-(n- ri)n(p) + iDq 2 + ir in ](p) 
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Substituting this form into Eq. (12) yields the following expression for the diffuson contribution to 
the correlator of LDOS 



Po(e) = 



T 



16tt 3 



l/v F Ttr 



dqq 

Jo 



2w 



S(k F \(f)p - 4>p>\) 



2tt Jo 2tt -i [e - (f, - rj)n(p)] + Dq 2 + r^(p) ' 



(28) 



Two integrations (over q and over p /) in Eq. (28) can be readily performed resulting in the following 
expression for SO-induced satellites Po(e) 



Po(e) 



where the functions C^ v (e) are the following azimuthal averages 



(29) 



o 2tt 



tr 



(30) 



Since C^^{e] 



£r),fj,{—£) the cooperon contribution, -Pc(e), to the LDOS correlator coincides with 
-Po(e), as it should be expected on general grounds. Thus, Eqs. (29) and (30) constitute our final 
result. The functions £ MjT/ (e), which determine the energy dependence of the correlator, exhibit a 
singular dependence on a weak magnetic field, oj z <C u so for /j, = —i], as demonstrated in the next 
subsection. 



D. Shape of the satellites 

As it follows from Eq. (26), in a weak magnetic field we have Tj nt (0) ~ 2r tr = const(cf)). Then for 
the satellites, centered at e = ±2uj so we can identify the ^-dependence of the integrand in Eq. (30). 
This dependence comes from fl(4>), defined by Eq. (5). For uo z <C u so , we have fl(4>) ~ uj so + uj z sin <fi. 
Upon substituting this form into Eq. (30), the integration can be performed analytically with the 
use of the following identity 



1 

2tt Jo 



2tt 



— / d(fi\n (a + 6 sin 0) z 



In 



2 In 161 



+ Va 2 - b 2 



for 
for 



\a\ > \b\ 
\a\ < \b\ 



(31) 



yielding 



P(e) = - 



In le - 2cu so | r tr 



irE F T tr J \ lnuj z r tr 



for 
for 



2^so| > ^z 

2a> so < 



(32) 



The satellites of the correlator -P(e), calculated numerically without the assumption u z <C uj so , are 
shown in Fig. 2. We see that the plateau in the shape of the satellites is very well pronounced, but 
slightly tilted. The reason for this tilt is the effect of the "background" originating from the main 
peak. 



E. Smooth potential; |ei — £2! <C uso 

In this subsection we solve the system of Eqs. (17) with kernels defined by Eq.(18)) for small 
energy difference |ei — e 2 | <C cu so . Firstly, we point out that the vertex functions and are 
relevant, since in Eq. (18) £1 must be equal to £ 2 - At the first glance, in equation (17) for T^, only 



G 



the term in the r.h.s. containing the same should be kept. Indeed, the coupling of this term 

to is determined by the kernel KZ^, which, as it is seen from Eq. (18) contains a square 

of the small parameter $- At)At . However, unlike the case of the satellites, the "feedback" becomes 
important for small |ei — e 2 |. Namely, the equation (17) for contains in the r.h.s. the coupling 

to Tfi£ with the same small coefficient oc Q 2 -^- Specifics of the small energy difference, as compared 
to satellites, is that both Ffr£ and TZfi£ are resonant in this case. Thus Eq. (17) reduces to the 
following system of coupled equations 

r£>, p>, q, U ) = S{k F \<P p - <f) pl \) + J^Kjfiip, Pl , q, u)r^( Pl ,p>, q, u ) 

+ j^M~4(p^^-ZM^^ ( 33 ) 



T-ZM,qM = S{k F \<f> P - p H)|<Wp,Pi)| 2 + J^KZ^{p^q^-Z^P'^M 

+ /^^b^i-^^r^b!,^^^). (34) 

Solution of the system (33), (34) yields [20] 

T^(p,p', q,uj)= ^Sihl^ - M) { 2 \. - 1( S + I n 2 lo- -if J > (3 5 ) 

w 2r [to + iDq 1 + ir f (p) tu + iDq 1 + 2ir ini (p) J 



uj + iDq 2 + irj 1 (p) 



u + iDq 2 + 2irJ(p) 



(36) 



where the diffusion pole in the first term is cut by the inelastic time, Tf. Substituting and r_ /(/( 



into the general expression Eq. 
at the final result 

P(e) = - 



[12) and keeping only the terms with coefficients $ w ~ 1, we arrive 



1 



2nE F T t r 




(37) 



Note that each term in the sum (37) is comprised of two (equal) contributions from T{ 1 and TZ^Zx- 
In the derivation of Eq. (37), as well as of Eq. (29), we used (f F r tr ) _1 as the upper cutoff in the 
integrals over q. In principle, there is also a contribution from the ballistic region, up to q ~ k F . 
This contribution does not affect the correlation function of the tunneling conductance, studied in 
Sec. IV. 



III. DISCUSSION 



The general case of an arbitrary SO-splitting would require a lengthy calculation, involving 4x4 
diffuson and cooperon matrices [18,19]. However, in the most interesting case of a strong splitting, 
when uj so r int ^> 1, a considerable simplification occurs: the LDOS correlation function, P(e), at large 
enough e can be expressed in terms of that in the absence of splitting, P(e,cu so — 0) = -Po(e) 

The splitting energy 2u so enters only through the argument of Pq and the low-energy cutoff h/T int , 
which assumes the role of dephasing time. 
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To establish the relation between P and Pq, we notice [20] that the states with the same energy, 
E, have two different values of the wave number, hk± = \j2mE ± lu so /v f , depending on the branch 
(see Fig. 1). Similarly, states at energy E' = E + e have wave numbers Tik'± = \/2mE + (e ± u so ) /v F . 
Thus, we can identify four contributions to P(w), coming from correlating each of the k± states 
with either of k'± states. Since, for u so r int ^> 1, the degree of correlation between a pair of states 
depends only on the difference of the wave numbers, all four contributions have the same functional 
form. It is given by the function Po( e ) with the corresponding arguments, which can be either 
Hv F (k' + — k + ) = Hv F (k'_ — fc_) = e or hv F (k' + — kJ) = e + 2cu so , or hv F (k'_ — k + ) = e — 2u so . It 
must be emphasized that, since separate branches cannot be resolved beyond the time T int , the 
correlator Pq for |e| ~ 2u so must be computed with a low-energy cutoff h/Ti n t. Hence the condition 
tTint ^> 1, adopted in this qualitative consideration, insures that the SO-induced satellites in P(e) 
do not overlap with the main peak. Combining the four contributions, we obtain 

P(e) = i [2P (e) + Po(e + 2co so ) + P (e - 2cu so )] . (38) 

Obviously, the above general consideration, applies to the both diffuson and cooperon contributions 
to the LDOS correlator, P. This consideration clarifies why the height of the central peak in P(e) 
is twice the height of each of the SO satellites. The factor 2 in the first term of Eq. (38) appears 
because the argument e occurs twice out of the four possibilities mentioned above. 



IV. MANIFESTATION IN THE TUNNELING SPECTROSCOPY 



In the tunneling-spectroscopy experiments [21,22,10,11,23,24,12-14] the measured quantity is a 
resonant-tunneling current, I sd , through a localized impurity state as a function of the source-drain 
bias, V^ d . Such an impurity plays the role of a "spectrometer" since its energy position, eo, changes 
with V sd . Within a narrow range, 5V, of V sd this change is linear, i.e., A = e (V + 5V) — e {V) ~ (38V, 
where (3 is the structure-specific parameter. Experimentally [21,22,10,11,23,24,12-14], measurements 
are performed in the plateau regime when the temperature is much lower than V^ d , so that the value 
of resonant-tunneling current 



I(v ) = £ r de lii^rje) 

ad ' hJ-oo [ e _ eo (v: d )p + [ 7 ,( e ) +7r ( e )] 2 /4 

is temperature independent. In Eq. (39), 7; and 7 r stand for the tunneling widths, associated 
with the escape from the impurity into the emitter and collector, respectively. These widths are 
proportional to the LDOS in the electrodes. Typically, the widths 7; and 7 r differ strongly, 7 r ^> 7;. 
In addition, the energy dependence of 7 r is weak, so that it can be considered as a constant. As 
a result, it is a tunneling coupling to collector that dominates the width of the Lorentzian in Eq. 
(39). With regard to the weak fluctuating dependence of current, /, on V sd in the plateau regime, it 
is exclusively due to the energy dependence of 7;, which, in turn, originates from the fluctuations of 
the LDOS. 

Quantitative analysis of the experimental data is performed (see e.g. Ref. [13]) by plotting the 
correlator C(5V) = (5g(V+5V)5g(V)) of the fluctuations of the differential conductance g = dI ad /dV 
around its average value, (g), as a function of 5V for different values of V = V sd . The expression for 
C(5V) in terms of the LDOS correlator P{e — e') is obtained (see [9,25]) by, first, calculating the the 
correlator of current fluctuations (81 (V sd ) 51 (V a ' d )) and then taking derivatives with respect to V sd 
and V' d . Using the fact that a convolution of two lorentzians is also a lorentzian, the final expression 
for C(5V) can be written in the form 
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C(SV) = (5g(V d + 5V)5g(V d )) 



W a 2 

4 <9A 2 



duo [P(c 



P(-uj) 



( w -A)' + 7= 



(40) 



If, following Ref. [9], we substitute the conventional form P(u) oc In (w 2 t 2 + l) into Eq. (40), where 
Tf is the "floating up" time of a hole created as a result of the tunneling act, then Eq. (40) yields 
[9] C(5V) = F 2 (5V/V C ), where the dimensionless function F 2 (x) is defined as 



Fo(x) 



l-x 2 
(1 + x 2 ) 



2 ' 



(41) 



and the characteristic value, V c , is given by V c = (3~ l (j r + fi/r t ). In the presence of the SO-coupling 
the correlator P(uo) has three peaks, centered at uo = and uo = ±2uo so . Without Zeeman splitting 
all three peaks have the same shape, which allows to express the correlator, Eq. (40), in terms of 
the function F 2 as follows 



/35V - 2co 3 

Wc 



+ -F 2 

4 2 



(35V + 2u s 



(42) 



where V' c = V c + P^h/n 



int 



f3-l( lr + h/T f + k/T mt ), V*' 



V c + 2(3~ l n/T, 



int 

It is 



(3~ x (7 r + h/r f + 2h/r int ), and r int is the intersubband scattering time defined by Eq. (26). 
seen from Eq. (42) that the SO peaks in P(oo) give rise to the satellites in C(5V) at 5V = 2f3~ 1 uo so . 
Thus, the SO-coupling has a noticeable effect on the correlator of the tunnel conductances when 
2(3~ 1 uj so > V' c . In Fig. 3 the correlator C is plotted vs. 5V jV' c for three different values of the 



dimensionless SO-splitting oj' sc 



,I(3V' C . It is seen, that at c^ G = 0.5 the correlator C (SV/Vf:) 



develops a characteristic "shoulder", while at lj' so = 0.7 it exhibits a well-developed additional maxi- 
mum. In the experiments [11-14] the characteristic width of the impurity level, 7 r , was < IK. This 
suggests that the SO satellites should be well resolved in tunneling spectroscopy of the narrow-gap 
semiconductors, where u; so is of the order of several meV. Concerning the Ga^s-based structures, the 
SO splitting there is much smaller (of the order of IK), and depends strongly on the details of the 
confinement potential. Therefore, the observation of the additional peaks in tunneling spectroscopy 
would require a higher "spectrometer resolution". The characteristic signature of the SO-related 
features in the correlator C(5V) would be a strong sensitivity of the correlator to a weak parallel 
magnetic field. In Fig. 4 we plot C (5V/V/,) for the same values u' so = 0.5 and uj' so = 0.7 as in Fig. 
3 and for different ratios u z /u so . The curves were calculated using Eqs. (29), (37) for P(cu). We see 
that both the "shoulder" at u' so = 0.5 and additional maximum at u>' so = 0.7 do not disappear with 
increasing u z , but rather get shifted down. Remarkably, the effect of the parallel field is noticeable 
already at very small u z /u so ~ 0.1. 



V. CONCLUSIONS 



In the present paper we have demonstrated that the intrinsic SO-splitting of the electron spectrum 
in the absence of disorder manifests itself in a disorder-induced effect - mesoscopic fluctuations of 
the local density of states. This observation suggests that the splitting, 2u so , must show up in 
the resonant tunneling spectroscopy in the form of additional peaks in the correlator C(5V). The 
fact that the positions of these peaks do not depend on the disorder offers a possibility to measure 
experimentally the magnitude of the splitting. Note that, unlike the magnetotransport experiments, 
where the SO-related features are quickly washed out with increasing temperature [26], the resonant 
tunneling current in the plateau regime depends on temperature rather weakly [12]. As it was 
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demonstrated in Sec. II, the correlator, P, of the LDOS develops a plateau in a weak parallel magnetic 
field u z <C u so . This singular behavior manifests itself in the anomalous sensitivity of the correlator 
of the tunneling conductances to u> z , as illustrated in Fig. 4. 
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FIG. 1. Schematic illustration of the processes responsible for SO-induced peaks in LDOS correlator at e = ±2w s 




FIG. 2. Correlator of the LDOS, P(ei,e2), calculated from Eqs. (29), (37) for uJsoT tr = 10 and Tf/r tr = 10, is plotted 
versus dimensionless energy e/u! S o- The numbers near the lines show the value of Uz/uso- For convenience different curves are 
shifted along the vertical axis. 
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FIG. 3. The correlator C(SV) = (8g(V)Sg(V + 8V)) of differential tunneling conductanaces is plotted from Eq. (42) versus 
dimensionless voltage difference SV/V^ for different values of the SO-splitting: u! SO = (solid line); lv so = 0.5V C ' (dotted line); 
uso = 0.7Vc (dashed line). For simplicity, we have assumed 7 r > h/n n t, so that V' c « V c ~ V c ". 




FIG. 4. Evolution of the SO-related features in the correlator C(SV) with parallel magnetic field is shown for same 
Uso = 0.5V C ' (a) and ut S o = 0.7V C ' (b) as in Fig. 3. Solid lines Uz = 0; dahsed lines lo z = O.lc^so; dotted lines co z — 0.2u! SO . 
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